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We generalize recent work on parametric resonances for nonlinear Schrodinger (NLS) type equa- 
tions to the case of three dimensional Bose-Einstein condensates at zero temperatures. We show the 
possibility of such resonances in the three-dimensional case, using a moment method and numerical 
simulations. 
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I. INTRODUCTION 



There has been a blossoming of literature on the fea- 
tures of systems exhibiting Bose-Einstein condensation 
(BEC), triggered by its recent experimental realization 
Intially experiments with 10 3 to 10 6 atoms of ru- 
bidium or sodium (later experiments have used lithium 
and eventually (spin-polarized) hydrogen), in harmonic 
or cigar shaped traps have demonstrated condensation to 
a "pseudo- macroscopic" level of occupancy of the ground 
state for nK temperatures. Time of flight measurements, 
velocity distributions as well as spatial profiles have con- 
vincingly supported the physical picture of an abrupt 
transition in the behavior of the Bose gas, which has 
been interpreted as the signature of BEC. 

Following these experiments, many theoretical studies 
were launched to characterize different aspects of Bose 
condensates such as hydrodynamic modes 0, collective 
excitations || , the behavior of ideal quantum fluids [ 4 , 5 1 , 
the fraction of noncondensate vs. condensate atoms [ 3, 7\, 
or the generation and stability of vortices ||. In turn, 
experimental studies have progressed to address some of 
the theoretical predictions ||] and open up new questions. 

Here, we concern ourselves with one aspect of these 
quantum fluids, namely parametric driving. For the pur- 
pose of this report, we will restrict ourselves to the frame- 
work of the mean-field or Hartree-Fock approximation. 
This approximation is rigorously justifiable only at T = 
but it is expected || that the contribution of the non- 
condensate to the density is quite small. It is well-known 
that at this mean-field level the condensate wavefunc- 
tion is governed by the Gross-Pitaevskii (GP) [jl0| equa- 
tion. An issue addressed after the original experiments 
achieving the condensation was the study of collective 
excitations ||(llj . In these papers these excitations were 
induced by a harmonic trap weakly modulated in time 
with appropriate types of symmetry. More recently, it 
was demonstrated that extended parametric resonances 



can occur |12| in a two-dimensional (2d) NLS equation 
with a harmonic trap. This result may or may not (for 
reasons to be explained below) be relevant for two di- 
mensional studies of Bose gases. However, this naturally 
raises the question of whether a similar result can be de- 



duced for the 3d case which is certainly of direct relevance 
to experimental studies. 

The main question we will address is whether weak har- 
monic modulation of trapped 3d condensates can cause 
an anomalously large response in their wavefunction. 
Our answer, which will be in the affirmative, will be moti- 
vated by mathematical analysis using a moment method 
and verified by numerical simulation. We will briefly dis- 
cuss the implications of these results and the suggestion 
of relevant experiments. 



II. MOMENT METHODS 

Considering a spherical trap, the dimensionless GP 
equation for the dynamics of the BEC condensate is Jl3] 
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Here, the subscript t denotes time derivative and £ = 
{%TtN\a\/ a^) 1 ^ is a dimensionless parameter arising 
from the number N of particles, the s-wave scattering 
length a and from ai characterizing the strength of the 
trap (see, Ref. p3|)- In Eq.(|l|) \(t) is a dimensionless 
function allowing for time dependence of the trap and 
v = sign(a), generalizes the equation to describe attrac- 
tive (a < 0) as well as repulsive (a > 0) interactions. 
Since, we restrict ourselves to spherical symmetry we 
only include the radial contribution in the Laplace op- 
erator. Although we are mainly interested in the full 
three-dimensional (3d) case we will in general consider 
the d-dimensional version of Eq.(§) so that 
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Similarly to Ref. JT^] , we define the following quantities 
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where (d) indexes the dimension. This type of nonlin- 
ear Schrodinger (NLS) equation (|l]) has two conserved 
quantities: as a result of the phase invariance the norm 
corresponding to j|1 q is conserved in any dimension d. 
Also, since Eq.(Q) is a Hamiltonian system arising from 



if = 



C 4 |Vm| 2 + -M 4 + A(i)r 2 |n| 2 J r^dr (7) 



this quantity is conserved. It is useful to note that the 
Hamiltonian (or more appropriately energy functional), 
if , can be expressed in terms of the moments Eqs.(^)-(|rj) 
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The relevance of the moments Eqs.(^)-(|^), is based on 
their time evolution and in the following we will derive 
the relations governing this dynamics. The physical ra- 
tionale behind such an approach lies in the fact that the 
resulting equations can yield predictive diagnostics for 
the dynamics of the BEC. In particular, will essen- 
tially yield the width of the spatial profile of the wave- 
function. If a condensation phenomenon is to take place 
even when starting from a spatially uniform distribution 
(at temperatures T < T c ), the width must evolve towards 
a constant non-zero value. 

Using Eq.(Q) and its complex conjugate one can derive 
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In deriving these we assume u to vanish as r — * oo. 

Unfortunately, it is in general impossible to close this 
hierarchy of equations because the time derivative cou- 



ples to the next order e.g. i. 
couples to 1^ 
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and so on. However, 



combining Eqs. (Q) and (|T^) provides some insight 
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First, this confirms that the norm I^q is conserved in any 
dimension. Secondly, this relation clearly suggests a = 2 
as good choice since the term involving ij a— 4 th en van- 
ishes irrespective of dimension. Also, this choice allows 
the use of Eq.rtq) to reduce the expression (113) to 



16C" 4 A(t)4 d 2 ) + 2v(d 
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which corresponds to the relation commonly referred to 
as the virial theorem |l4[ for the nonlinear Schodinger 
without a trap, A(£) = 0. 

Clearly, the 2d case is special as the structure of 
Eq.(|l2"|) is such that a closed time evolution can be pre- 
scribed. In addition, for time dependent modulation of 
the trap amplitude, we find (as was observed for this 
problem in Ref. |lj| and again in Ref. jL2j) a Hill type 
equation which establishes parametric resonances for the 
behavior of the width (or amplitude) of the wavefunc- 
tion. Important as this conclusion about the 2d behavior 
may be for general NLS-GP equations, it is not clear 
that it is relevant to BEC. Since BEC is not possible 
in spatial dimensions less that three [d < 3) || where 
a Kosterlitz-Thouless topological transition seems to be 
occurring instead |5|]l6fl , the applicability of the GP equa- 
tion for d < 3 is controversial. Our search for condensate 
instabilities is therefore most compelling in three dimen- 
sions where no such reservations exist. 



Although Eq.(]12j) is not closed in 3d, it is easily seen 
that the following inequalities hold for d > 2 
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where the equality applies to the two-dimensional case 
only (and in fact for the noninteracting 3<i case v = 0). 
The value of these inequalities lies in the predictions 
about the 3d case ]1^| . Since we can resolve, or at least 
very well characterize, the 2d behavior, we are now able 
to extend this to quantitative predictions about of the 
3d behavior. For the attractive case v < the possi- 
bility of collapse occurs as /| 2 can become zero in fi- 
nite time. In 2d and due to Eq.([l4"|) also in 3d, a suf- 
ficient condition for collapse is H < 0, although, de- 
pending on the initial configurations, collapse can be 
achieved even for if > 0. A more complete discus- 
sion is given in Ref. frL8| . For the more realistic case 
(in BEC contexts) of repulsive interaction v > we see 
for example that the parametric resonances that were 
demonstrated in the two-dimensional case Q will also 
be present in the three-dimensional case. So, for instance, 
for d = 3 and X(t) = Ag(l + ecos(wt)) will exhibit res- 
onance around Ao = nto/2, n = 1,2,3.. where the ex- 
tent of the first resonance is determined by the inequality 
|1 — ix) 2 /(4Ag)| < e/2 However, additional resonances 
may be possible in 3c? due to the nonlinear driving result- 
ing from the repulsive interaction. Some understanding 
of the influence of the last term in Eq. fll2] ) on the dynam- 
ics can be gained by assuming that the wavefunction can 
be approximated as 



2 



5(K 




small. In the repulsive v > case however B will gen- 
erally not become small since there is no collapse. This 
simple analysis suggests that the parametric forcing of 
the experimentally realizable 3d case will result in a res- 
onance picture analogous to that previously reported for 
the 2d problem. Our numerical simulations of the full 
Eq.(|l]) with A(t) as given above verifies the validity of 
this prediction. A typical example for Aq = 1, u> = 1, 
(i.e. n = 2) and e = 0.05 is given in Fig.l. The re- 
sponse of the wavefunction (whose initial condition had 
max. x \tji(x, t = 0)| 2 = 1), corresponding to the parametric 
resonance, can be observed directly from the wave func- 
tion Fig. 1(a) but even more clearly in the time evolution 
of I2 2 , as shown in Fig. 1(b). 



III. CONCLUSION AND FUTURE CHALLENGES 
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1. (a) evolution (actually r \u(r, t)\ | is plotted for 
clarity) of the condensate wavefunction both in space (r) and 
in time (i) for part of the domain (close to r = 0) and (b) 
the evolution if of i^a- Parameters are Ao = 1, w = 1, (i.e. 
second resonance) and e — 0.05. 



r(3) 



-B- 3 / 2 ip(r/B), 
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where U\ = r 2 ip(r)dr is a shape-dependent constant. 
Thus, assuming adiabatically the wavefunction does not 
alter its shape ip significantly as a result of the dynamics, 
u as defined in Eq. (U) automatically satisfies the norm 
conservation. Utilizing this in Eq.(|l2|) yields 



d 2 B 2 
dt 2 



\{t)B 2 =Qi + vQ 2 B- 



(16) 



where Q\ and Q2 are constants determined by the shape 
if) and the initial value of B. Clearly, the last term in 
Eq. ( |l6| ) will only influence the dynamics when B becomes 



In this paper, we have presented and extended the for- 
malism of the moment method, used in Refs. |0| and |2(J 
for the 2d GP equations, to the more relevant 3d case. 
We have commented on the special nature of the two- 
dimensional problem where the moment equations form 
a closed set of equations. We have also added a note of 
caution in considering the results of the GP analysis for 
d < 3. It might well be that, analogous to mean-field 
analysis in statistical physics systems, the "critical di- 
mension" for this system is, indeed, d c — 3 and for lower 
dimensionalities the predictions of the mean-field theory 
are unreliable. A satisfactory self-consistent first prin- 
ciples description of an interacting boson gas for d < 3 
presents a very challenging theoretical problem, and it re- 
mains an unresolved issue whether a transition is present 
(and if it is, what is its nature). 

On the other hand, we have used the moment methods 
and have derived results for GP functionals in all dimen- 
sions of physical interest. Considering, in particular, the 
3d case, where the validity of the GP approximation is 
clear, we have obtained a non-closed set of equations for 
the moments of the wavefunction. We have demonstrated 
that for a parametric time-disturbance of the trap am- 
plitude, parametric resonances are possible. To date the 
experiments that have used parametric modulation have 
not observed such phenomena. These experiments have 
been performed in cigar-shaped traps (where the analysis 
is considerably more complicated even in the 2d problem 
]l2t ) . No fine tuning of frequencies and amplitudes was 
explored since the aim of the studies was to excite collec- 
tive modes rather than to observe parametric resonances. 
Hence, we propose an experiment in a spherical trap us- 
ing weak harmonic modulation of the condensate. Given 
the current experimental advances (see e.g. Ref. J2IJ for 
a review), such an experiment seems feasible. Such a 
study would, apart from the validation of the theoreti- 
cal prediction, also explore how such resonances might 
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destabilize the condensate. 
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